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ABSTRACT 

The purpose of this letter is to check the quality of different methods for estimating stellar masses of galaxies. 
We compare the results of (a) fitting stellar population synthesis models to broad band colors from SDSS and 
2MASS, (b) the analysis of spectroscopic features of SDSS galaxies (Kauffmann et al. 2003), and, (c) a simple 
dynamical mass estimate based on SDSS velocity dispersions and effective radii. Knowing that all three methods 
can have significant biases, a comparison can help to establish their (relative) reliability. In this way, one can also 
probe the quality of the observationally cheap broadband color mass estimators for galaxies at higher redshift. 
Generally, masses based on broad-band colors and spectroscopic features agree reasonably well, with a rms scatter 
of only ~ 0.2 dex over almost 4 decades in mass. However, as may be expected, systematic differences do 
exist and have an amplitude of ~ 0.15 dex, corrleting with Ha emission strength. Interestingly, masses from 
broad-band color fitting are in better agreement with dynamical masses than masses based on the analysis of 
spectroscopic features. In addition, the differences between the latter and the dynamical masses correlate with 
Ha equivalent width, while this much less the case for the broad-band masses. We conclude that broad band 
color mass estimators, provided they are based on a large enough wavelength coverage and use an appropriate 
range of ages, metallicities and dust extinctions, can yield fairly reliable stellar masses for galaxies. This is a very 
encouraging result as such mass estimates are very likely the only ones available at significant redshifts for some 
time to come. 

Subject headings: galaxies: mass function — galaxies: fundamental parameters 



1. INTRODUCTION 

The stellar mass of galaxies at the present epoch and the 
build-up of stellar mass over cosmic time has become the fo- 
cus of intense research in the past few years. 

In the local universe, results on the stellar mass function of 
galaxies were published using the new generation of wide-angle 
surveys in the optical (Sloan Digital Sky Survey; SDSS, York 
et al. 2000; 2dF, e.g. Folkes et al. 1999) and near-infrared (Two 
Micron All Sky Survey; 2MASS, Skrutskie et al. 1997). Cole 
et al. (2001) combined data from 2MASS and 2dF to derive the 
local stellar mass function, Bell et al. (2003) used the SDSS and 
2MASS to the same end. 

At z > 0, a number of authors studied the stellar mass density 
as a function of redshift (Brinchmann & Ellis 2000; Drory et al. 
2001; Cohen 2002; Dickinson et al. 2003; Fontana et al. 2003; 
Rudnick et al. 2003) reaching z ~ 3, while others, using wider 
field surveys, investigated the evolution of the mass function of 
galaxies (Drory et al. 2004; Fontana et al. 2004) to z ~ 1.5. 

Generally, the high-redshift work relies on fits of multi-color 
photometry to a grid of composite stellar population (CSP) 
models to determine a stellar mass-to-light ratio, since large and 
complete spectroscopic samples of galaxies are not yet avail- 
able. A similar approach was chosen by Cole et al. (2001) and 
Bell et al. (2003), too, at z ~ 0. 

Taking advantage of the availability of photometry and spec- 
troscopy for galaxies in the SDSS, Kauffmann et al. (2003, K03 
hereafter) utilized spectroscopic diagnostics (4000A Break, 
D n 4000, and the H5a Balmer absorption line index ) to estimate 
the mean stellar age and the fraction of stars formed in recent 



bursts in each galaxy. By comparison of the colors predicted 
by their best-fit model to the object's broad-band photometry 
they determine the amount of extinction by dust and hence the 
stellar mass-to-light ratio. 

The purpose of this letter is to compare the stellar masses 
determined by this spectroscopic technique to masses obtained 
from multi-passband photometry and to compare both methods 
to a simple dynamical estimate of mass. Knowing that none 
of these methods yields a fiducial (stellar) mass, a comparison 
helps to establish the (relative) reliability of each method and 
makes us aware of potential differences between these estima- 
tors. Moreover, it can show us whether one can use observa- 
tionally cheaper estimators as surrogates for more expensive 
(or unobtainable) ones, which is particularly important when 
dealing with high-redshift datasets. 

Specifically, we want to know how the two estimators com- 
pare to each other, if using K-band M/L yields better masses 
than using the g-band M/L which is accessible at high z, and 
how these compare to a simple dynamical mass estimator, 
M ~ <7 2 R e /G. 

This letter is laid out as follows. In Sect. 2 we describe the 
sample of galaxies we use in this work. In Sect. 3 we give a 
brief overview of how we derive stellar masses by fitting CSP 
models to multi-band photometry. In Sect. 4 we compare these 
masses to the values in K03 and in Sect. 5 to a simple dynam- 
ical estimate of mass based on the SDSS velocity dispersions. 
We also discuss the implications of these comparisons. 

We assume £2 M = 0.3, Q. A = 0.7, and H = 70 km s -1 Mpc -1 
throughout this work. 
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FIG. 1. — Illustration of the model-fitting technique used to estimate stellar masses. The lower panels show the comparison of the best fitting model to the 
photometric data. The red and blue lines represent the main and burst component, respectively. The green line represents the combined SED. The left hand side 
shows a young object with a burst component, the right side an older object. The upper panels show projections of the likelihood function onto four planes, age vs. 
dust in the main component, burst fraction vs. burst extinction, age of the main component vs. burst fraction, and star formation timescale, T, vs. burst fraction. The 
resulting likelihood distributions of M/L are shown in the upper right panel on each side (M/L g blue; M/L,- green; M/L/c red). The M/L of the best fitting model is 
indicated by vertical lines. 



2. THE GALAXY SAMPLE 

The sample of galaxies we use in this work is selected from 
the NYU Value-Added Galaxy Catalog 1 (Blanton et al. 2004). 
This is a merged catalog of objects from the SDSS Data Release 
Two (DR2) and 2MASS point-source and extended-source cat- 
alogs (and other catalogs which are not relevant here, as well). 

We select all objects classified as galaxies and having a se- 
cure redshift measurement in the SDSS and that are detected in 
the 2MASS catalogs. From this set we randomly sub-select 
20% of the objects leaving us with a sample of sample of 
~ 17000 objects having redshifts and photometry in ugrizJHK. 

We cross-correlate this catalog with the data from K03 2 to 
obtain their stellar mass estimates. 

The galaxies in the sample span the absolute magnitude 
range —15.3 < M g < —23.5, the restframe u — g color range 
0.5 < u—g < 2.0, and the (stellar) mass range 8 < logM < 12. 

3. DERIVING STELLAR MASSES 

The method we use to infer stellar masses from multi-color 
photometry is an advancement of the program used in Drory 
et al. (2004). It is based on the comparison of multi-color pho- 
tometry to a grid of stellar population synthesis models cover- 
ing a wide range in parameters, especially star formation histo- 
ries (SFHs). 

We base our new model grid on the Bruzual & Chariot (2003) 
stellar population synthesis package. We parameterize the pos- 
sible SFHs by a two-component model, consisting of a main 
component with a smooth analytically described SFH and a 
burst of star formation. The main component is parameterized 
by a star formation rate of the form y/(f) exp(— t/t), with 
T G [0.1, oo] Gyr and a metaUicity of -0.6 < [Fe/H] < 0.3. The 

1 see also http : / / was sup .physics .nyu. edu/vagc/ 

2 available online at http : / /www . mpa-garching . mpg . de/SDSS/ 



age, f, is allowed to vary between 0.5 Gyr and the age of the 
universe (at the object's redshift). 

The smooth component is linearly combined with a burst of 
star formation, which is modeled as a 100 Myr old constant star 
formation rate episode of solar metallicity. We restrict the burst 
fraction, j3, to the range < j3 < 0.15 in mass (higher values 
of j3 are degenerate and unnecessary since this case is covered 
by models with a young main component). We adopt a Salpeter 
initial mass function for both components, with lower and up- 
per mass cutoffs of 0.1 and 100 M Q . 

Additionally, both the main component and the burst are al- 
lowed to exhibit a variable amount of extinction by dust. This 
takes into account the fact that young stars are found in dusty 
environments and that the starlight from the galaxy as a whole 
may be reddened by a (geometry dependent) different amount. 
In fact, Stasihska et al. (2004) find that the extinction derived 
from the Balmer decrement in the SDSS sample is independent 
of inclination, which, on the other hand, is driving global ex- 
tinction (see, e.g., Tully et al. 1998). This is different from the 
approach taken by K03, where a single extinction value for the 
whole galaxy is used. 

We compute the full likelihood distribution on a grid in 
this 6-dimensional parameter space (r, [Fe/H],f ,Ay,f},Ay), the 
likelihood of each model being °< exp(— ;£ 2 /2). To compute 
the likelihood distribution of M/L, we weight the M/L of each 
model by its likelihood and marginalize over all parameters. 
The uncertainty in M/L is obtained from the width of this dis- 
tribution. 

This procedure is illustrated in Fig. 1, where we show SEDs 
and likelihood functions for two objects, a young object with a 
high burst fraction, and an older and fairly quiescent object. We 
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FIG. 2. — Comparison of our photometry-based stellar mass estimates to the stellar masses of K03. The left panel shows the masses of K03 plotted against our 
estimate based on M/L g , the right hand panel against masses based on M/Lk- The colors denote Ha equivalent width from no emission (red) to strong emission 
(> 25A; purple). The small panels show histograms of the residuals again as a function of Ha equivalent width. 



show projections of the likelihood function onto four planes in 
parameter space, age vs. dust in the main component, burst frac- 
tion vs. burst extinction, age of the main component vs. burst 
fraction, and star formation timescale, t, vs. burst fraction. The 
figure also shows the resulting likelihood distributions of M/L 
in the g, i, and K bands. Note that for the quiescent object, 
the width of the M/L distribution is very similar in the g and K 
bands, while it is much wider in g than it is in K for the younger 
star forming object. On average, the width of the likelihood dis- 
tribution of M/L at 68% confidence level is between ±0.1 and 
±0.2 dex (using M/Lg). The uncertainty in mass has a weak 
dependence on mass (increasing with lower S/N photometry) 
and much of the variation is in spectral type: early-type galax- 
ies have more tightly constrained masses than late types (see 
also Fig. 1). Using the U band, the uncertainty in mass grows 
by -0.05 dex. 

4. COMPARISON OF STELLAR MASS ESTIMATORS 

In Fig. 2 we compare our photometry -based stellar mass esti- 
mates to the stellar masses of K03. We show the Ha equivalent 
width (as measured by the SDSS) by color coding. The overall 
impression from this figure is that the two different estimators 
agree remarkably well, within a rms scatter of only — 0.2 dex 
over almost 4 decades in mass (and hence they largely agree 
within their respective uncertainties). This is only a relative 
statement, though. It does not imply that the masses are ac- 
curate to that level in an absolute sense, although it is very 
reaffirming. However, as may be expected, there are system- 
atic differences as a function of star formation activity on the 
— 0.15 dex level. 

The D n 4000 and H5a based method of K03 yields masses 
almost identical to ours for weakly star forming objects 
(EqW(Ha) < 5A) at masses above 10 M©. At lower masses, 
our estimator tends to give slightly higher masses than K03's. 
For more strongly star forming objects, the photometrically de- 
termined masses are smaller than the ones of K03. For objects 
with EqW(Ha) > 25A, the discrepancy becomes as large as 
0.15 dex, independent of mass. Note that at high redshift, such 
objects will be more common. 



We suspect that there are a multitude of reasons for these 
differences based on the different sampling of stellar popula- 
tions by both methods (if we leave out the JHK bands, our 
masses become more similar in their trends to K03's, although 
with increased scatter). Plausibly, though, this is explainable by 
the fact that the photometrically determined masses sample the 
light from the whole galaxy, while K03's SDSS-based sampling 
of D n 4000 and H5a covers only the inner 3 arcseconds. Since 
most galaxies are redder in their centers than in the outer parts, 
this might lead to higher masses for star forming disk galaxies. 
Early-type galaxies without blue star forming disks do not suf- 
fer from this effect. This is confirmed by restricting the sample 
to low redshifts, which maximizes the effect. Also, at the lowest 
masses, galaxies might have more irregular SFHs, and photo- 
metric methods might fail in this case (Bell & de Jong 2001). 
However, Fig. 2 does not show a dependence of the residuals on 
mass, only on current star formation rate. 

Fig. 2 also shows that our masses based on the g band are 
very similar to the ones estimated through the K band. For early 
type systems and weakly star forming systems, they are statisti- 
cally indistinguishable. Star forming systems, however, tend to 
have g-band masses lower by — 0. 1 dex. The good agreement 
is partly due to the fact that the effect of dust extinction and 
age on the effective M/L are very similar in canonical models, 
as has been pointed out by Bell & de Jong (2001), although 
they are not completely degenerate. Typically, the spread of the 
stellar M/L in the galaxy population at any given luminosity is 
around 0.7 dex in g and 0.35 dex in K. These results is reaf- 
firming since the restframe blue spectral range is accessible to 
photometry to very high redshift, and thus high-z studies mostly 
rely on M/L B . 

5. COMPARISON WITH DYNAMICAL MASSES 

Since we do not have a fiducial mass estimator (neither for 
stellar mass nor for total mass, for that matter), it is only nat- 
ural to ask how the stellar mass measurements presented here 
compare with estimators of mass based on kinematic data. In 
fact, it has been suggested that stellar mass (or, more accurately, 
baryonic mass) and total mass are tightly related and that stel- 
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FIG. 3. — Comparison of stellar mass vs. dynamical mass, Md yn °= a 2 R e /G. The left hand panel shows M^ yn vs. photometric mass and the right hand panel shows 
Mjy,, vs. the masses of K03. Color encodes Ha equivalent width, as in Fig. 2. 



lar mass can be used as a surrogate for total mass in the context 
of high-z galaxy surveys to probe structure formation (see, e.g., 
Brinchmann & Ellis 2000). 

We use the measurements of velocity dispersion, a, and ef- 
fective radius in the g band, R e , provided by the SDSS pipeline 
to plot stellar mass vs. dynamical mass, Md yn « a 2 R e /G, in 
Fig. 3. The left hand panel shows Md yn vs. photometric stellar 
mass and the right hand panel shows Mdyn vs. the stellar masses 
of K03. Color again encodes Ha equivalent width. 

Above 10 10 M Q , the stellar masses from both methods fol- 
low the dynamical masses remarkably well. Below ~ 10 9 Mq, 
the velocity dispersion measurements of the SDSS becomes un- 
reliable as we approach the instrumental resolution of the data 
(-TOkms- 1 ). 

At higher masses, although both estimators generally follow 
Mdy n , there are again some differences, and both estimators 
show similar trends in their residuals although with different 
amplitudes. Stellar masses agree very well with Md yn at the 
highest masses (which are mostly populated by old, quiescent 
objects). At lower masses, stellar masses show a trend to larger 
values than Md yn with decreasing mass and with increasing Ha 
emission line equivalent width. This effect is weak in the pho- 
tometric estimator, and stronger in K03's method, which gives 
stellar masses larger than Md yn by 0. 1 to 0.4 dex at almost all 
masses. This comparison is unchanged by using M/Lk instead 
ofM/L g . 

It is important to note that Md yn is not a good estimator of 
total mass, and that this comparison is again only to be taken in 
relative terms. In fact a 2 R e /G can only provide a lower limit to 
the mass. However, as long as a bulge is present, the total mass 



should not be underestimated by more than ~ 0.3 dex (see, e.g. 
Fig. 4 in Whitmore & Kirshner 1981; also Padmanabhan et al. 
2004, who show that a 2 R e /G is a reasonable mass estimator, 
although this paper is concerned with ellipticals only). It is 
therefore not surprising to find stellar masses in excess of Md yn 
and the difference between the two increasing at lower masses. 

Nevertheless, the point of this work is to assess the general 
consistency and reliability of stellar mass estimates than to in- 
vestigate the relationship between stellar and dynamical mass 
in galaxies. Fig. 3 shows that the estimators of stellar mass, and 
especially the photometric estimator which is most easily ob- 
tainable for large high redshift samples (covering a bluer wave- 
length range, though), closely follow Md yn as measured by this 
simple dynamical measure. We cannot see significant system- 
atic deviations which would bias or invalidate this estimator. 
This is a very encouraging result, since such an estimator is 
very likely to be the only one available at z > for some time 
to come. 
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